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Abstract 

We use the technique of conformal transformations to generate self-similar collapse 
in Brans-Dicke theory. We analyze the solutions concerning the critical behavior found 
recently by Choptuik. The critical exponent associated to the formation of black hole 
for near critical evolution is obtained. The role of the coupling parameter is discussed. 

PACS number: xxx 



Recently, Choptuikjl] has showed numerically the occurrence of critical phenomena in 
spherically symmetric gravitational collapse of massless scalar field. These critical phenom- 
ena are characterized by power-law behavior, scaling relations and a form of universality. 
In his study, Choptuik analyzed the evolution one-parameter families of initial conditions 
(parameter p, say), and classified the solutions in three classes: (i)subcritical (p « p cr ), 
where the final state is a complete dispersion of scalar field resulting in flat spacetime; (ii) 
supercritical (p >> p cr ) in which the scalar field is intense enough to trap itself in a black 
hole; (iii)critical solutions (p = p cr ), that represents the transition between complete disper- 
sal and black hole formation. The most interesting results merge for near critical (p = p cr ) 
evolutions due to two important features. First, the mass of formed black hole in this regime 
is given by M BH oc \p — p cr | 7 , where 7 = 0.37 is conjectured to be universal since it does not 
depend on the particular initial condition. Second, for near critical evolution a strong field 
region exists where the scalar field oscillates in a very particular way. Such oscillations are 
echoes and the scalar field satisfies 0(e~ n A t, e~ n A r) = </>(£, r), where is an integer number 
and A = 3.4. It is important to mention that the exact critical solution has an infinite train 
of echoes in the strong field region. 

Quite surprisingly the same results were obtained for axisymmetric collapse of gravita- 
tional waves and for the spherically symmetric collapse of radiative fluids || exhibiting 
local self-similarity. Therefore, it is strongly suggested that these critical phenomena are 
independent of the symmetries as well as the collapsing matter. Analytical solutions de- 
scribing such features could be extremely useful for their understanding. Unfortunately, this 
task has proved to be very difficult, and an alternative approach was to suppose continuous 
self-similarity. In several self-similar exact solutions with ordinary scalar field the critical ex- 
ponent was found to be 0.5||], whereas for the conformally coupled case, 0.21||. Maison||, 
analyzing the collapse of perfect fluid, showed the non-universality of the critical exponent 
in the sense that it depends strongly on the state equation relating the energy density and 
pressure. 

In this paper, our objective is to study analytically the self-similar collapse in Brans- 
Dicke theory of gravitation (BDT), and its relation with the critical behavior as described 
above. We are motivated by the renewed interest in the called scalar-tensor theories of 
gravitation in which the BDT is the most simple. We mention, for instance, the numerical 
work performed by Scheel et al|| concerning the collapse of matter in BDT. We will consider 
here a simplified model where the matter field is absent. Then, we expect to find, if it is 
present, the role of the coupling parameter u on the value of the critical exponent. 

BDT is the most simple of the scalar-tensor theories. In such theories the spacetime is 
characterized by the metric tensor g^ v and by a scalar field coupled to the metric and matter. 
The action for the BDT is: 

S = J d^x^g-^R-^g^d^d^ + Lj (1) 

where is the scalar field, u the coupling constant and L m is the matter Lagrangian. For 
a more general scalar-tensor theory u = to{(p). The basic characteristic of this theory is 
the variable gravitational term measured locally as G = ^ ■ According with the most 
recent experimental tests, u has minimum value about 500. As showed by Dicke||, it is 
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possible to express the theory in the Einstein frame. In this frame the action (1) is formally 
given as General Relativity with ordinary scalar field, which is a suitable redefinition of the 
scalar field 6. Then, we have: 



S 



(fxJ-g (R - g^ u <9 M $ d v $) 



where R is the scalar of curvature obtained from the metric g^ u , $ is the ordinary scalar 
field and we have considered L m = 0. Both actions integrals are equivalent if the following 
conformal transformation holds: 



1 „ 

and the scalar field $ is expressed in terms of (f> by: 



l2u + 3 



In (ft 



(3) 



(4) 



Therefore, from the above equations, it is possible to generate solutions in BDT starting 
from those solutions for ordinary scalar fields in GR. 

To generate solutions in BDT describing self-similar spacetimes we consider the Roberts' [10 
solution. Such a solution has been studied as an analytical model that exhibits critical be- 
havior in scalar field collapse M. The Roberts' solution is given by: 



ds 2 = —d udv + r 2 (u, v) (d 9 2 + sin 2 9 d <p 2 ) 
where u and v are null coordinates, and 

1 



r (u,v) 



— — + a 2 u + b 2 v 
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with: 



f± 



4 a 2 u — (1 ± \ 1 — 16 a 2 b 



v. 



(5) 



(6) 



(7) 



The scalar field $ is: 



±—= In 
V2 



f+(u,v) 



The collapse situation is properly described if a 2 is chosen positive, more specifically, a 2 = 
1/4, without loss of generality. The solutions for which a 2 < corresponds to expanding self- 
similar cosmological models. The collapse solutions are classified according to three distinct 
values of the parameter b 2 : < b 2 < 1/4, the subcritical regime, b 2 = 0, the critical regime 
and the supercritical b 2 < 0. In this last case, there is the formation of a "black hole" type 
configuration (Fig. 1(a)), meaning that the formed singularity is enclosed by an apparent 
horizon, but the spacetime is not asymptotically flat. The critical exponent related to near 
critical evolution was found to be equal to 0.5. 

To obtain the spherically symmetric self-similar collapse in BDT, we begin with Eqs. (4) 
and (8), which yields: 
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where = ^2^+3 • According with Eqs. (3) and (9), we have: 

ds 2 = -(p^dudv + T^iu.v) {dd 2 + sin 2 6d V 2 ) (10) 
where the proper area, £ 2 («, v), is given by: 

S 2 (n, v) = J- i;)- 1 /_(«, u) /+(«, v). (11) 

ID Q2 

The conformal factor, 0, is regular everywhere with exception of f-(u, v) = or f+(u, v) = 0, 
depending of the sign of (3(u). Indeed these regions represent, in general, singularities^. As 
an important property of conformal transformations, the causal structure of the spacetime 
is not changed if the transformation is well behaved, which is the present case. We as- 
sume uj > 0, implying |/3(o>)| < 1 and, as a consequence, the singular regions have zero 
proper area. Otherwise, for uj < (|/3(<^)| > 1), the singular regions could be character- 
ized by an infinity proper area. Two types of solution emerge from Eq. (9): the one with 
positive sign inside the parenthesis, and another with negative sign. For the first type, 
we must have ^ + > 0. Hence, we chose a 2 = 1/4 so that the proper area is given by 

£(«, v ) 2 = \ f-(u, v) 1+l3 ^ f + (u, d) 1- ^^. For the second type, it is necessary that ^ + |^| < 0. 
Consequently, a 2 must be negative in order the proper area to be positive as required. How- 
ever, due to the fact the collapse is not described for such a choice of a 2 , we are not going to 
study this situation. Nonetheless, a 2 can be positive only if |/3(u;)| is an even number that 
is realized only by uo < 0. 

The apparent horizon is determined locally by: 

<T £ „£„ = (). (12) 
After substituting Eq. (11) into (12), two relations are obtained: 



u AH = (l + (3(uj)Jl-4b 2 )v (13) 



46 2 f , iJlN 

U AH = ~ , > r j-^ (14) 

where the first equation corresponds to £ )U = whereas the second £ „ = 0. Note that in 
the limit uj —> 00, we recover the results provided by Roberts' solution as expected. The 
next step is to verify for which cases the apparent horizon is present. 

As mentioned, the Penrose diagrams are not altered by the conformal transformation. 
Therefore, following ref. ||, we classify the solutions as subcritical (0 < b 2 < 1/4), critical 

/ ± jl-2/fi \ 

1 For instance, the scalar of curvature is given by R = 16 u> 1 (1—4 62) u V \ j ■ Since w > 0, \j3\ < 1 

meaning that at f—(u,v) — or f + (u,v) = 0, the scalar of curvature diverges. For sake of completeness, 
only if \(3\ > 2 (or -1.5 < uj < -1.375) R is finite at f-(u,v) = (fl < 0) and f+(u,v) =0(J3> 0). 
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(p2 = 0) and supercritical (62 < 0). Studying the presence of the apparent horizon in the 
collapse in BDT the following results are obtained. For £> 2 = 0, it is not difficult to show 
that there is no apparent horizon described either by Eq. (13) or Eq. (14), unless uj < — 1 
that has no interest. The subcritical case (0 < 6 2 < 1/4) displays some novelty. Taking 
(3{uj) = ^2ij+3 < 0' an d after a straight analysis, we conclude that there is no apparent 
horizon. However, for (3{uj) = ^ 2 ^ +3 > 0, the apparent horizon given by Eq. (14) exists 
for any value of to > and encloses the timelike singularity. This situation is also found for 
the collapse of conformally coupled scalar field ||, but has no relevance regarding the critical 
behavior. 

The most interesting case is the supercritical (62 < 0). Let us consider first f3(u) > 0, 
where, after easy verification, the apparent horizon indicated by the line OH (Fig. 1(b)) is 
spacelike and described if by Eq. (14) with the condition uj > uj* = — 2 6 2 — 1 holds. For 
uj — > 00 we recover the apparent horizon given by Roberts' solution (line OR, Fig. 1(b)). On 
the other hand if uj = uj*, the horizon becomes the null surface v — 0. In Fig. 2, we depict 
the effect of varying uj on the apparent horizon. Finally, considering /3(uj) < the apparent 
horizon is again described by Eq. (14), but with no restriction on possible values of uj. The 
effect of varying uj from —3/2 to 00 is to rotate the line OH from the f-axis to the line OR 
(the apparent horizon of Roberts' solution) as indicated in Fig. 2. 

Black Hole mass and the Power Law for Near Critical Behavior 

We are going to see the influence of the coupling parameter in the power law associated 
to the mass of the formed black hole. This last quantity is the effective gravitational mass 
inside the apparent horizon: 

m AH = 1 y/f_(u,vy+K»)f+(u,vy-ft") \ah (15) 

Substituting Eq. (14) with 62 < into the above equation and considering near critical 
evolution, 6 2 — 0, we have: 

m AH 3 (1 - g( "' ) ~ 4M [-(3 - m) t, + 1 - fa)]*** (-^ v. (16) 
1 — p{u>) 

As mentioned previously, self-similar spacetimes fails to be asymptotically flat and black 
hole configurations have infinite mass when v — > 00. However, an important point to be 
stressed is the obtained power-law txiah — (— ^) 7 for near critical evolution. Notice that the 
exponent 7 depends strongly on uj, or, in another word, of the intensity of which the scalar 
field couples to the geometry. Then, as showed in Fig. 3, if (3{uj) > 0, we have 0.79 < 7 < 0.5 
for uj varying from to 00. In the case of (3{uj) < 0, 0.21 < 7 < 0.5 for uj varying from to 
00. Incidentally, the exact value 0.37 is achieved in this case when uj = 5.8. A large value 
of the coupling constant compatible with earlier experimental results, uj = 500, produces 
7 = 0.52 and 0.48 for f3{yj) positive and negative, respectively. 

We can argue that the above result as well as the one obtained for conformally coupled 
scalar field rise doubts regarding the universality of the value 0.37. In these models, however, 
there a serious weakness which is the assumption of continuous self-similarity, and contrary 
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to the problem treated by Choptuik, the spacetime is not asymptotically flat. Koike and 
Mishimaflll]] studying the collapse of a thin massive shell coupled with an outgoing null 
fluid resulting in an asymptotically flat model showed that the critical exponent is not 
universal. On the other hand, from the large amount of works dealing with self-similar 
collapse in connection to cosmic censorship hypothesis, a cutoff is introduced to add an 
asymptotically flat region. This procedure is equivalent to truncate the spacetime at some, 
say, v = Vq. Hence, we can match the self-similar spacetime with Schwarzschild or outgoing 
Vaidya solutions, where in the later case it models the radiation that escapes during the 
collapse. In general, a null shell characterized by mass and surface pressure [PH would be 
present for such a junction . In this way, the local trapped mass, represented by the mass 
function at the apparent horizon, is given by the mass of the shell and the Brans-Dicke scalar 
field mass. A reasonable condition to be imposed is that the shell has negligible mass and 
pressure, which could be seem as natural from the physical point of view. If this hypothesis 
is satisfied, Eq. (16) is rewritten as: 

1 — U\U))) 2 1-/3M 1+PM 

m AH - 1 P J H 3 - ^ 62 + 1 - * H*)~^ «b- ( 17 ) 

1 — p\u) 

We expect that, depending on the asymptotically flat spacetime to be joined, the conditions 
m s heii = and P s heii — at the apparent horizon, render some relation between vq and 
hi- Then, it will be possible to estimate, in this simplified model, the value of the critical 
exponent. We are investigating these issues and the results will be given elsewhere ||13||. 
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• Fig. 1(a) "Black hole" type configuration. The spacetime is not asymptotically flat 
and the mass of the black hole grows without bound as v — > oo. At v = we match 
smoothly the scalar field solution with the Minkowski spacetime. The lines indicated by 
R and H represent the apparent horizon of the Roberts and the Brans-Dicke solutions, 
repectively. 

• Fig. 1(b) The lines OR and OH represent the apparent horizon for the Brans-Dicke 
(case (3{uS) > 0) and Roberts solutions, respectively. The effect of varying u from cj* 
to oo is to rotate OH in the counterclockwise sense. For u = oo both lines coincide. 

• Fig. 2 For the case f3(uj) < 0, the apparent horizon behaves in a different manner as 
compared to the previous case. As far as u grows, the line OH rotates in the clockwise 
sense and tends to OR for u — > oo. 

• Fig. 3 Behavior of the critical exponent 7(0;) = ^±|H . 
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